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SUMMARY 


This  report  describes  modeling  methods  that  allow  the  computation  of  point  predictions 
and  prediction  probability  intervals  for  cumulative  workers’  compensation  costs.  Underlying 
these  models  is  the  actuarial  loss  development  factor  method,  a  method  that  computes  projected 
costs  by  utilizing  ratios  of  known  cumulative  costs  in  consecutive  years.  While  the  relationship 
between  cumulative  loss  development,  cohort,  and  development  year  in  these  models  is 
nonlinear,  a  transformation  renders  them  in  the  form  of  standard  linear  statistical  models,  thus 
allowing  the  development  of  prediction  probability  intervals  when  the  error  structure  is 
Gaussian.  The  modeling  methods  are  illustrated  using  data  collected  from  U.S.  Department  of 
the  Navy  workers’  compensation  payments  made  from  1990  through  1993,  including  claim  costs 
originating  from  1961  through  1993. 
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1.  Introduction 


The  usual  and  preferred  actuarial  method  for  projecting  workers’  compensation  claim 
costs  is  based  on  the  so-called  Loss  Development  Model.  To  illustrate  this  method.  Table  1.1 
lists  actual  total  (indemnity  and  medical  combined)  cumulative  claim  costs  (in  thousands  of 
dollars)  from  cohorts  of  claimants  originating  in  years  1990  through  1993.^ 

The  column  for  development  year  1  is  simply  the  total  claim  costs  incurred  for  the  year 
in  which  the  claims  originated,  namely,  the  cohort  year.  Cumulative  costs  are  available  for  the 
1990  cohort  for  years  1990,  1991,  1992,  and  1993.  In  contrast,  data  are  not  available,  for 
example,  for  the  1992  cohort  in  development  years  3  and  4,  since  these  would  be  accumulated  in 
the  years  1994  and  1995  and  are  not  yet  available. 


Table  1.1.  Actual  Total  (Indemnity  and  Medical)  Cumulative  Claim  Costs  (In  Thousands  $) 

Development  Year 


Cohort  Year 

i 

2 

3 

4 

1990 

14,955 

41,424 

62,897 

79,971 

1991 

13,566 

40,314 

60,137 

♦ 

1992 

14,468 

41,892 

* 

* 

1993 

13,702 

♦ 

* 

* 

It  is  of  interest,  based  on  the  above  data,  to  project  (or  forecast,  predict)  the  missing  costs 
represented  by  the  asterisks.  Actuaries  approach  this  problem  by  computing  Loss  Development 
Factors  (LDF)  for  consecutive  years.  For  a  given  cohort,  these  are  simply  the  ratios  of  the 
cumulative  costs  from  consecutive  years.  Actuaries  treat  LDFs  across  cohorts  as  homogeneous 
and  average  them  with  the  intention  of  improving  precision.  Finally,  they  compute  cumulative 
LDFs  by  taking  cumulative  products  of  these  averages.  These  computations  are  summarized  for 
the  data  in  Table  1.2.  Notice  that  by  construction,  the  successive  costs  for  a  given  cohort  may  be 
computed  by  multiplying  the  cohort  year  cost  (development  year  1)  successively  by  the  LDFs. 
For  example,  the  1990  cohort  generates  (in  thousands)  $14,955  the  first  year,  ($14,955) 
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(2.76991)  =  $41,424  accumulated  through  the  second  year,  ($14,955)  (2.76991)  (1.51837)  = 
$62,897  accumulated  through  the  third  year,  and  so  on. 


Table  1.2.  Loss  Development  Factors  Computed  from  Table  1.1. 


Cohort  Year 

Develooment  Year 

1  to  2  2  to  3 

3  to  4 

1990 

2.76991 

1.51837 

1.27146 

1991 

2.97169 

1.49172 

- 

1992 

2.89549 

- 

- 

Average: 

2.87903 

1.50504 

1.27146 

Cumulative: 

2.87903 

4.33307 

5.50932 

Through  experience,  actuaries  have  found  that  applying  the  estimated  (average)  LDFs 
calculated  in  this  fashion  to  the  cohorts  with  missing  data  yields  accurate  projections.  That  is,  to 
project  the  accumulated  cost  through  development  year  2  for  the  1993  cohort,  one  simply 
computes  ($13,702)  (2.87903)  =  $39,449.  Similarly,  the  accumulated  cost  (using  cumulative 
LDFs  calculated  from  the  average  LDFs)  through  development  year  3  for  the  1993  cohort  is 
projected  to  be  ($13,702)  (4.33307)  =  $59,372.  Proceeding  in  this  fashion,  the  missing  values  in 
Table  1.1  are  projected  as  shown  in  Table  1.3. 


Table  1.3.  Point  Predictions  of  Missing  Costs  from  Table  1.1. 


Development  Year 

Cohort  Year 

i 

2 

3 

4 

1990 

- 

* 

- 

- 

1991 

- 

- 

- 

76,462 

1992 

- 

- 

63,049 

80,165 

1993 

- 

39,449 

59,372 

75,489 

Actuaries  use  numerous  variations  of  this  method.  For  instance,  the  method  may  be 
applied  separately  to  indemnity  and  medical  costs  to  achieve  greater  projection  accuracy.  Also, 

! 

I 

j 
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volume-weighted  averages  may  be  used  instead  of  simple  averages  of  the  columns  of  LDFs. 
Finally,  the  costs  entered  in  the  model  may  or  may  not  be  adjusted  for  economic  factors,  such  as 
inflation,  with  the  understanding  that  the  resulting  predictions  would  retain  the  same 
interpretation  with  respect  to  these  factors.  That  is,  if  the  costs  entered  into  the  model  are  (or  are 
not)  adjusted  for  inflation,  then  the  predictions  are  (or  are  not)  adjusted  for  inflation.  Typically, 
however,  actuaries  only  present  their  results  as  point  predictions,  and  no  method  for  computing 
measures  of  prediction  accuracy  are  offered.  (Indeed,  searches  of  the  literature  and  consultation 
with  practicing  actuaries  indicate  that  such  methods  have  not  been  widely  investigated.) 

Many  sources  of  potential  error  are  inherent  in  the  LDF  method.  Data  of  this  nature 
contain  random  fluctuations  and  other  errors  (e.g.,  systematic,  specification).  Also,  when  a  large 
number  of  cohort  years  are  involved,  comparable  LDFs  across  cohorts  can  have  a  trend.  For 
example,  the  ratio  of  development  year  2  cumulative  cost  to  development  year  1  cost  for  the 
1961  cohort  could  be  significantly  larger  than  that  for  the  1992  cohort.  Having  not  accounted 
for  this  trend,  its  effects  show  up  as  error  in  the  basic  LDF  model  previously  described.  To 
develop  an  assessment  of  prediction  accuracy,  the  LDF  method  must  be  generalized  and  put  into 
the  context  of  a  statistical  model  that  accounts  for  and  makes  assumptions  concerning  the  nature 
of  statistical  error.  The  natural  way  in  which  to  assess  prediction  error  is  to  accept  the  prediction 
itself  as  a  random  variable  having  a  certain  probability  distribution,  and  to  provide  not  only  a 
numerical  value  for  the  prediction,  but  also  an  interval  with  the  interpretation  that  the  true  value 
of  the  cost  will  fall  into  that  interval  with  a  preselected  probability. 

Mathematically,  the  problem  may  be  described  as  follows.  A  cost  prediction  C  is  a 
function  of  available  observed  data  Cj,  C2, say  C  =  f(Cj, ...,  C^^).  A  prediction  interval 

(L(Cp  ...,  C^,  U(Cj,  ...,  Cjj))  for  C  with  coverage  probability  p  is  a  random  interval  with  the 

property  that  p  =  P|  L(Cj, ...,  C^^)  <  C  <  U(Cj, ...,  C^^)  }.  Naturally,  the  larger  (or  smaller)  the 

level  of  p,  the  wider  (or  narrower)  the  prediction  interval  will  be.  For  a  fixed  p,  the  width  of  the 
prediction  interval  reflects  the  accuracy  of  the  prediction. 

The  remainder  of  this  report  develops  intrinsically  linear  statistical  LDF  models  from 
which  predictions  and  prediction  intervals  may  be  computed.  The  basic  models  are  nonlinear. 
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but  are  termed  "intrinsically  linear"  because  they  become  standard  linear  statistical  models  after 
a  natural-log  transformation.  This  intrinsically  linear  structure  makes  the  development  of 
prediction  probability  intervals  more  tractable.  Section  2  develops  the  mathematical  structure  of 
the  models,  and  section  3  applies  the  models  to  the  computation  of  predictions  and  prediction 
intervals  for  U.S.  Department  of  the  Navy  workers’  compensation  costs. 

Two  sources  of  data  are  used  in  section  3.  Complete  data  collected  by  the  U.S. 
Department  of  the  Navy  for  1990,  1991,  1992,  and  1993  cohorts  are  used  to  illustrate  one  of  the 
models.  Yearly  incremental  costs  and  claim  counts  were  not  available  in  the  U.S.  Department  of 
the  Navy  database  for  the  1961  through  1988  cohorts  for  development  years  prior  to  1990  (but 
are  available  from  1990  forward).  To  study  the  available  data  and  to  make  future  year 

3 

projections  of  total  cumulative  costs  Miccolis  used  actuarial  methods  to  "reverse  forecast"  the 
cumulative  costs  for  the  1961  through  1989  cohorts  for  development  years  prior  to  1990  by 
examining  trends  in  the  incremental  cost  and  claim  count  data  available.  These  imputed  data, 
along  with  the  actual  data,  are  used  by  Miccolis  to  then  produce  projections  of  future  year 
cumulative  costs.  In  section  3,  these  imputed  and  actual  data  are  used  as  though  they  were  all 
actual  data  to  illustrate  the  use  of  one  of  the  models  developed  herein. 

The  final  section  contains  conclusions  of  this  investigation  and  recommendations  for 
further  study. 

2.  Intrinsically  Linear  Statistical  Loss  Development  Models 

Throughout  this  section,  the  basic  data  available  will  be  assumed  to  consist  of  cumulative 

costs  (accumulated  over  development  years)  C-  -,  1  <  i  <  N-j+1,  1  <  j  <  N.  Here,  the  subscript  i 

*0 

designates  the  ith  cohort  and  the  subscript  j  designates  the  development  year.  N  is  the  total 
number  of  cohorts.  Usually,  and  this  will  be  assumed  here,  the  cohorts  are  from  successive 
years  and  arranged  in  increasing  order  by  year,  and  assigned  values  1  through  N,  respectively. 
This  facilitates  analysis  of  a  trend  in  LDF  across  cohorts.  In  the  example  in  the  introduction,  N 
=  4,  and  the  cohorts  1990,  1991,  1992,  1993  are  assigned  values  1,  2,  3,  and  4,  respectively. 
Only  the  upper  left  triangular  (including  the  main  diagonal-see  Table  1.1)  area  of  the  cost 
matrix  elements  are  observed;  the  object  is  to  predict  the  cost  values  in  the  lower  right  triangle 
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(and  beyond). 

Define  Y..  =  log(C.  ./C-  •  ,),  for  1  <  i  <  N-j+1,  2  <  j  <  N.  Notice  that  exp(Y.;)  is  the 

LDF  relating  year  j  cumulative  costs  to  year  j-1  cumulative  costs  for  cohort  i.  It  is  postulated 
that 


=  f(i  j)  +  £: 


y 


(2.1) 


where  f(i,j)  is  an  appropriate  deterministic  function  of  the  cohort  year  index  and  the 
development  year  j.  For  the  time  being,  the  error  terms  e-  are  assumed  to  be  independent 

2 

Gaussian  random  variables  with  common  zero  mean  and  finite  variance  o  Wj  >  0,  allowed  to 
depend  on  j  in  such  a  way  that 


i; 


,OQ 

.  -W.<oo. 

0=1  J 


(2.2) 


In  (2.2),  it  is  assumed  that  the  weights  Wj  >  0  are  known,  but  the  parameter  >  0  is  unknown. 


The  exact  form  of  the  function  f  in  (2.1)  is  unknown,  but  approximations  are  easily 

derived  based  on  plausible  analysis  and  empirical  evidence.  For  a  fixed  cohort  i,  the  LDFs  must 

approach  1  as  development  year  j  increases.  This  is  true  because  eventually  the  claims  for  each 

claimant  in  the  ith  cohort  must  cease.  Factors  affecting  this  include  mortality,  claim  settlement, 

and  injury  recovery.  Empirical  evidence,  such  as  data  collected  by  the  National  Council  on 

Compensation  Insurance,  suggests  that  for  fixed  i,  exp[f(i,J)]  should  decrease  monotonically 

and  smoothly  to  1  as  j-^.  Equivalently,  f(i,j)  should  decrease  monotonically  and  smoothly  to  0 

as  J-^oo.  To  illustrate  this.  Figure  2.1  shows  plots  of  LDFs  versus  year  for  accident  cohorts  from 

2 

1979  through  1990  as  reported  by  the  National  Council  on  Compensation  Insurance  from  data 
reported  by  private  insurers  providing  workers’  compensation  coverage  in  37  states.  These  LDF 
plots  are  for  indemnity,  medical,  and  total  costs.  Each  data  point  is  computed  from  a  five-year 
average. 
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Ist/and  2nd/3rd  3rd/4th  4th/5th  5th/6th  6thmh  7th/8th 

Year 
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Figure  2.1.  Nationwide  Loss  Development  Factors 
From  the  National  Council  on  Compensation  Insurance2 


Therefore,  regardless  of  the  exact  form  of  f.  it  can  be  postulated  that  for  fixed  i,  f  has  an 
asymptotic  expansion  (as  of  the  form 


f(id)  -  Aq  +  A/j  +  B/j^  +  C/j^  +  D/j"^  + . . . 


(2.3) 


(Here,  the  notation  g(j)  -  h(j)  as  means  that  lim.  g(j)/h(j)  =  1.)  The  leading  two  terms  in 

(2.3)  can  be  eliminated  as  follows.  Ultimately,  there  must  be  an  upper  bound  for  C-  •  as  j 

increases.  That  Is,  cohort  i  eventually  ceases  to  generate  further  costs.  The  model  (2.1)  entails, 
by  iteration,  that 


‘^i.j  =  ^i,l  nic=2“P<'^ik>  =  Ci,iexp(l[^2Y  J  =  exp(5;[^|f(i.k)  +  (2.4) 


Because  of  (2.2),  the  random  series  converges  (to  a  random  variable  that  is  finite)  with 

probability!.  By (2.3), 

XUl«M)  -  li=l[Ao  +  Aflc  +  B/k^  +  C/k^  +  D/k'*  +  . . .  ] 

and  since  Aq  j-x«  and  -  log(j)  -►«>  as  while  <  oo  for  m  >  1,  then 

limj^^i._jf(i,j)  will  be  finite  if  and  only  if  Aq  =  A  =  0.  Thus,  setting  Aq  =  A  =  0  in  (2.3)  and 

truncating  the  expansion  after  a  finite  number  of  terms  should  present  a  fairly  accurate 

approximation  of  f(i,j)  for  fixed  i.  To  verify  this  hypothesis.  Figure  2.1  shows  the  model  (2.3) 

2 

fit  to  the  data  from  the  National  Council  on  Compensation  Insurance  with  Aq  =  .A  =  0  and 
terms  beyond  order  4  neglected. 

From  Figure  2.1  it  may  be  concluded  that  qualitatively,  the  model  (2.3)  is  a  plausible 

2 

model  for  LDF  data.  It  is  noted  that  the  National  Council  on  Compensation  Insurance  data  is 
based  on  commercially  managed  payment  systems  that  are  very  different  from  the  mechanisms 
present  in  the  U.S.  Department  of  the  Navy.  In  particular,  the  rate  of  claim  settlement  is  much 
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higher  commercially.  Therefore,  it  would  not  necessarily  be  true  that  the  parameters  in  model 
(2.3)  that  apply  to  the  National  Council  on  Compensation  Insurance  data  would  be 
approximately  the  san»e  as  those  that  apply  to  the  U.S.  Department  of  the  Navy  data.  In  fact, 
with  much  slower  claim  settlement,  a  better  empirical  fit  of  (2.3)  tc  the  U.S.  Department  of  the 
Navy  data  may  be  achieved  by  including  the  terms  Aq  and  A/j,  and  relaxing  the  condition  (2.2). 

nile  theoretically  it  should  be  true  that  Aq  =  A  =  0  and  (2.2)  holds,  real  data  in  which  claim 

settlement  is  slow  (i.e.,  the  approach  of  f(i,j)  to  0  as  is  apparently  slow)  may  not  be 
extensive  enough  to  effect  a  good  fit  to  (2.3)  with  these  conditions  forced. 

The  study  of  the  variation  of  f(i,j)  for  fixed  j  is  based  mostly  on  empirical  observation. 
Evidence  from  both  Miccolis  and  the  National  Council  on  Compensation  Insurance  suggests 
there  is  a  slight  decreasing  trend  in  LDF  as  cohort  year  increases  for  fixed  j.  Figure  2.2  shows 
LDFs  from  Miccolis  for  j=2  varying  from  the  1961  cohort  to  the  1992  cohort.  A  quadratic 
trend  fit  by  least  squares  is  also  shown  in  Figure  2.2.  (It  is  noted  that  the  data  representing  the 
1961  through  1989  cohorts  were  imputed  in  Miccolis,  as  discussed  in  the  introduction.)  Of 
course,  as  j  increases,  the  total  variability  of  f(t,j)  diminishes,  so  there  is  the  need  to  allow 
interaction  (i.e.,  deviation  from  purely  additive  effects  of  the  variables  i  and  j)  between  the 
independent  variables  i  and  j. 
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Fitted  *  Data 


Figure  2.2.  1  st/2nd  LDF  Versus 
Cohort  From  1961  to  1962 


Assuming  f  is  a  smooth  function  of  i  and  1/j,  then  regardless  of  the  true  form  of  f,  a 
bivariate  Taylor  expansion  of  f  should  provide  an  approximation  (in  this  case  a  sum  of 
monomials  in  the  variables  i  and  1/j)  that  is  reasonably  accurate.  Retaining  terms  only  up  to 
order  4  in  1/j  and  order  2  in  i,  this  leads  to  the  consideration  of  a  parametric  model  for  f  of  the 
form 


1^=0  A 


(2.5) 


In  fitting  the  model  (2.5)  to  actual  data,  it  will  often  turn  out  that  not  all  terms  are  needed.  For 
example,  a  stepwise  regression  procedure  will  systematically  add  only  the  terms  necessary  in 
(2.5)  to  effect  a  good  fit  to  the  data  without  overfitting  or  underfitting.  For  further  discussion  of 

4 

over-  and  underfiiting,  see  the  report  by  Angus. 


2.1  Matrix  Formulation  of  the  Model 

The  LDF  model  developed  thus  far  may  be  expressed  as  a  standard  linear  model  as 
follows: 

First,  for  j  >  2  and  i  >  1,  let  v(i,j)^  denote  the  row  vector 


v(i,j)^=(  1  i  1/j  1/j^  1/j^ 


,,.4  ..-2  •r3  -A  -2,.  .2,.2 

1/j  i/j  i/j  i/j  i/j  1  /j  1  /j 


and  let  P  denote  the  column  vector 


3  =  ( 


"^10  ^20  "^01  "^02  "^03  ^04  ^11  ^12  ^\3  "^21  ^22  ^23  ^24 


)‘ 


with  signifying  the  transpose  of  the  matrix  M.  Also,  let 


’^'=(^12^22  "•  ^N-1,2 


^13  ^23  -  %-2,3  ••• 


^1,N-1  \n-1 


1,4 
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^~ri2^22  "  %-l,2 


^2,N-1 


and  X  be  the  n  by  15  matrix  whose  transpose  is  given  by 


X^  =  (v(l,2)  v(2,2) ...  v(N-l,2) 


v(l,3)  v(2,3) ...  v(N-2,3) 


v(l,N-l)  v(2,N-l) 


v(l 


,N) ), 


where 


n=N(N-l)/2 

is  the  total  number  of  components  of  Y.  Finally,  let  W  be  the  n  by  n  diagonal  matrix 

W  =  diag(w2  W2  ...  W2IW3  W3  ...  W3|...|wj^_j  '^N-iI'^n) 

where  the  first  block  containing  W2  is  of  length  N-1,  the  second  block  containing  W3  is  of  length 
N-2,  and  so  on.  Then  using  (2.5)  the  model  (2.1)  may  be  written  as 

Yij  =  v(i,j)'p  +  ejj  (2.6) 

or  in  matrix  notation  as 

Y  =  Xp  +  e.  (2.7) 

The  error  vector  in  (2.7)  is  transformed  to  one  with  homogeneous  variances  by  premultiplying 
by  to  yield 

W'^^^Y  =  W'^^^XP  +  e*  (2.8) 

2 

where  now  e*  is  n-variate  Gaussian  with  zero  mean  and  variance  matrix  o  I  (I  being  the  n  by  n 
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identity  matrix).  Standard  linear  statistical  model  theory,  such  as  that  found  in  the  text  by 
Arnold,^  now  provides  the  result  that  the  maximum  likelihood  and  best  linear  unbiased  estimator 
of  P  is  the  weighted  least  squares  estimator  given  by 


P  =  (XV^X)‘^XVW 


(2.9) 


and  an  unbiased  estimator  of  a  is 


y‘|w^-w^x(xV^x)‘^x*w^Jy 

=^(Y-Xpyw'^(Y-Xp) 


(2.10) 


where  p  is  the  number  of  unknown  parameters  in  p.  (It  is  possible  that  not  all  15  of  the  terms  in 
(2.5)  will  be  needed.  For  the  terms  not  needed,  the  corresponding  columns  of  X  and  parameters 
in  P  are  removed,  and  the  dimensions  of  X  and  P  adjusted  accordingly.  i  is  assumed  that  these 
adjustments  have  been  made  throughout  the  analysis.) 

Additional  facts  from  linear  statistical  model  theory  that  will  be  useful  in  developing 
prediction  intervals  are  that 

$  i  N(P,  o^(xV^X)'^),  (2.1  i) 

(n-p)oV  =  xVp),  (2.12) 

A  A2 

and  that  P  and  a  are  statistically  independent  Here,  N()i,  Z)  signifies  a  random  variable  that  is 

2 

normally  distributed  with  mean  (scalar  or  vector)  )i  and  variance  (scalar  or  matrix)  Z,  and  %  (m) 

signiHes  a  random  variable  with  a  chi-squared  distribution  with  m  degrees  of  freedom.  The 
d 

notation  X  =  Y  means  that  the  random  variables  X  and  Y  have  the  same  distribution. 

2.2  Point  Prediction 

The  first  step  in  prediction  of  costs  is  the  prediction  of  the  missing  log-LDFs.  Define 
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the  predicted  value  of  Y-j  by 


^  f  Y  .  if  Y  ;  is  observed; 

Y„  =  {  'J  (2.13a) 

[  v(i,j)^  P  otherwise  (j  >  2), 


so  that  the  predicted  LDF  going  form  year  j-1  to  year  j  for  cohort  i  is 
LDF. .  =  exp(Yjj) 


(2.13b) 


and  the  predicted  cumulative  LDF  (CLDF)  going  from  year  1  to  year  j  is 

CLDFij  =  lt.2^Fik  =  “PdL^ik)  (2.13c) 

(This  usage  of  the  term  "cumulative  LDF"  is  slightly  different  from  common  usage  in  actuarial 

2 

literature.  Typically,  as  defined  by  the  National  Council  on  Compensation  Insurance, 
cumulative  LDF  is  computed  from  a  given  base  year  to  an  "ultimate."  This  "ultimate" 
corresponds  to  the  limit  as  in  this  model.  This  approach  was  not  adopted  here,  since  no 
"ultimate"  costs  were  available  to  aid  in  fitting  the  models,  and  predictions  out  to  32  years,  as 
well  as  intermediate  years,  were  of  interest.) 

In  (2.13a),  P  is  the  estimator  given  in  (2.9).  Recalling  that  the  model  (2.6)  or  (2.7) 
entails  that 


(2.14) 


for  2  <  j  <  N-i+1,  and  treating  the  initial  costs  C.  ,  as  constants  (1  <  i  <  N),  a  natural  predictor 

for  missing  C- ;  is  given  by 
*0 


=  C^exp! 


(2.15) 


-15- 


with  Y. ,  given  by  (2.13a).  The  predictions  (2.15)  are  nearly  the  same  as  those  that  would  be 

lyX 

computed  using  the  actuarial  method  described  in  the  introduction  except  that  the  unknown  (lOg) 
LDFs  are  predicted  using  the  regression  equation  structure  (2.6)  rather  than  simple  averaging 
over  (nonhomogeneous)  cohorts.  Thus,  (2.15)  should  yield  more  accurate  predictions. 

To  illustrate  the  development  of  the  predictions  in  (2.15),  the  following  assumt  N  = 
4.  Here,  the  upper  left  triangular  (including  the  diagonal  elements)  numbers  ar  ually 
observed  and  do  not  need  to  be  predicted,  but  nevertheless,  (2.15)  reduces  to  the  actual  observed 
values.  The  lower  right  part  of  the  matrix  illustrates  exactly  which  (log)  LDFs  are  predicted 
using  the  estimated  regression  equation. 

^1,1  Cj  jexp(Yj2)  ^1,1®’^P^^12'^^13'^^14^ 

^2,1  ^2,1®^P^^22^  ^2,1®^P^^22‘‘‘^23^  ^2,1^^P^^22'*'^23''’^24^ 

^3,1  ^3,1®’^P^^32^  ^3,1®^P^^32'^^33^  ^3,1®^P^^32'^^33'*'^34^ 

^4,1  ^4,1®’^P^''''42^  ^4,1®^P^42'^^43^  ^4,1®^P^^42'^''43'^^44^ 

It  is  apparent  from  (2.13)  through  (2.15)  that  the  predicted  costs  also  satisfy  the  recursion 

CiJ  =  Cij_iexp(Yjj)  (2.16) 

where  j  ^  2,  which  is  exactly  the  "loss  development  principle;"  the  predicted  cumulative  cost  for 
year]  is  the  predicted  cumulative  cost  for  year  j-1  times  the  estimated  LDF  for  year  j-1  to  j. 

2.3  Interval  Prediction 

Suppose  that  in  addition  to  predicting  C-  •,  it  is  necessary  to  find  random  variables  L  and 
U  so  that  for  a  fixed,  preselected  probability  1-a,  P{L  <  C-  •  <  U}  =  1-a.  Here,  it  is  assumed  that 
i  ^  N  and  j  >  N-i+1.  (If  j  ^  N-i+1  then  C-  •  is  observed  and  there  is  nothing  to  predict)  The 
quantity  given  by 


(\  I/O 

is  pivotal;  that  is,  it  has  a  distribution  that  does  not  depend  on  any  unknown  parameters.  In 
(2.17),  o^  is  given  by  (2.10)  and  is  given  by  (2.13a). 

To  see  why  (2.17)  is  pivotal,  notice  that  Yy^,  k  >  N-i+2,  are  not  observed  in  the  current 
data  and  are  considered  future  observations.  Thus,  by  the  assumptions  underlying  (2.6),  it 

A 

follows  that  Yy^,  k  >  N-i+2,  are  independent  of  Yy^,  since  the  latter  are  functions  of  the  current 
data.  Since 


d  2  ^ 

where  Ey^  =  N(0,  o  Wj^)  is  independent  of  p,  it  follows  that 


and  hence  from  (2. 1 1) 


^k=N-i+2 


(Yik-V=N 


(2.18) 


where 


?  w. 

^k=N-i+2  k 


+ 


(xV*x)’^ 


(2.19) 


It  follows  from  independence  and  (2.12)  that 
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(2.20) 


that  is,  has  a  t-distribution  with  n-p  degrees  of  freedom.  Denoting  by  t^^Cm)  the  “U  quantile  of  the 
t-distribution  with  m  degrees  of  freedom,  (2.20)  implies  that 


A 


(2.21) 


which  in  turn  implies  that 


p{C;_jexp(T,.o  tj^2(ii-p))  <  Cj  j  <  Cj  jexp(t..o  t,.^(n-p))  |  =  1-a  (2.22) 


Relation  (2.22)  says  that  with  probability  1-a,  the  future  cost  C.  •  will  lie  in  the  interval 

*0 


^Cyexp(t;jO  t^2<"-P))  •  Ci  j«p(v  ‘l-o/2<“-rt))- 


(2.23a) 


In  (2.23a),  it  is  tacitly  assumed  that  0  <  a  <  1.  Notice  that  the  left  endpoint  is  the  point 

A 

prediction  C  .  multiplied  by  a  factor  less  than  1  (i.e.  reduced,  since  t  (m)  <  0  for  -u  <  1/2)  and 

A 

the  right  endpoint  is  the  point  prediction  C-  -  multiplied  by  a  factor  greater  than  1  (i.e.  increased, 

^0 


since  t^(ni)  >  0  for  “U  >  1/2). 


By  taking  C  ,  s  1  for  all  i,  the  analysis  that  led  to  expression  (2.33a)  also  provides  (1-a) 

A 

prediction  intervals  for  cumulative  LDF  via 


®^P&k=2Yik  ^ij^  'a/2^"-P0  ’  ^^PfLl^ik  ^ij^  h-a/2^"-P>)]-  ^2.23b) 

Similarly,  prediction  limits  for  individual  LDFs  are  easily  developed.  Letting 
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r?j  =  w.  +  v(i,j)‘(xV^X)'^(i.j) 


(2.23c) 


it  follows  that  a  1-a  prediction  interval  for  the  LDF  exp(Yjj)  is  (assuming  j  >  N-i+l) 

^exp(Yj.  +  rj.o  .  exp(Y.j  +  (2.23d) 

2.4  Geometric  Average  Model 

In  this  section  it  is  assumed  that  the  function  f(i,j)  in  (2.1)  does  not  depend  on  i  (i.e,  that 
LDFs  are  homogeneous  across  cohorts  but  may  change  with  development  year  j).  Thus,  assume 
that  for  all  i, 

f(ij)  =  HC).  (2.24) 

This  would  partly  justify  the  actuarial  method  of  prediction  discussed  in  the  introduction,  but 
instead  of  taking  simple  averages  of  like-LDFs  as  advocated  in  the  actuarial  method,  the 
following  will  make  a  case  for  using  geometric  averages  instead. 

Instead  of  parametric  modeling  of  the  variation  in  expected  log-LDF  over  j  as  in  the 
previous  sections,  it  is  also  possible  to  treat  each  value  of  n(»)  in  (2.24)  as  a  parameter  to  be 
estimated.  The  disadvantage  here  is  that  predictions  will  not  be  available  for  C.  for  j  >  N,  since 

observed  LDFs  will  not  be  available  in  this  range.  It  could  be  argued,  however,  that  prediction 
beyond  the  development  year  range  of  the  data  is  not  advisable  anyway. 

Another  difficulty  with  the  approach  of  this  section  is  that  there  is  loss  in  efficiency  by 

introducing  so  many  more  parameters.  In  the  previous  method,  no  more  than  15  parameters  (not 
2 

including  a  )  are  estimated  from  the  data  while  in  the  present  model,  there  must  be  N-1 
parameters  estimated  (namely  |i(2),  ...,  ii(N)).  Loosely  speaking,  for  N  >  16  and  to  achieve  a 
given  level  of  precision,  more  data  will  be  needed  for  the  model  ot  this  section  than  for  the 
regression  model  of  section  2.2. 
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Under  this  set  up,  the  parameters  p(j)  are  estimated  by 


A  .  1  Y'N-j+l 

"  N-j+l^r=l 


A  ♦ 


Predictions  C-  •  of  unknown  costs  C-  •,  2  <  j  <  N,  are  constructed  as  before  using 
i»J  ij 


except  the  predicted  Y-jS  are  now  computed  via 


A  ♦ 


Yjj  if  Yjj  is  observed; 

^l(j)  otherwise  (2  <  j  <  N). 


As  before,  a  prediction  of  LDF  from  year  j-1  to  year  j  for  cohort  i  becomes 


A  ♦  A  ♦ 

LDF.j  =  exp(Y.j) 


and  the  predicted  cumulative  LDF  (CLDF)  going  from  year  1  to  year  j  is 


The  parameter  is  now  unbiasedly  estimated  by 


—  2  2  y^N-1  1  Y^N-j+1  /  .  \2 

^  ^  (N-l)(N-2)  ^j=2  ^  ^i=l  (^ij  ’  ■ 


The  quantity 


(2.25) 


(2.26) 


(2.27a) 


(2.27b) 


(2.27c) 


(2.28) 
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A  * 


(Vik 


(2.29) 


is  pivotal.  To  see  this,  standard  Gaussian  linear  model  theory  applies  allowing  one  to  conclude 
that 


yj 

Mc=N-i+2 


A  ♦ 


(2.30) 


where 


2  _  yj  w  + 

“  '^k=N-i+2'*'k  ik=N-i 


w, 


N-i+2  N-k+1 


(2.31) 


and  that 


(N-l)(N-2)5^d  2 
2 - Ly  J..  =  X  ((N-1)(N-2)/2) 

2a^ 


independently  of  (2.30),  so  that 


(2.32) 


(2.33) 


where  n  =  N(N-l)/2  as  before. 

At  this  point,  it  is  noteworthy  to  compare  (2.33)  with  (2.20).  In  (2.20)  the  degrees  of 
freedom  are  n-p  with  p  <  15,  regardless  of  the  size  of  N,  while  in  (2.33)  the  degrees  of  freedom 
are  n-(N-l),  which  eventually  (as  N  grows)  is  much  smaller  than  n-p.  Thus,  (2.33)  will  be  more 
variable  than  (2.20)  and  lead  to  wider  prediction  intervals  than  (2.20).  On  the  other  hand,  for  N 
<  6,  n-p  can  be  less  than  or  equal  to  zero  and  then  the  parametric  method  based  on  (2.5)  would 
not  be  applicable  (that  is,  there  are  too  many  parameters  and  not  enough  data).  In  the  present 
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model,  n-(N-l)  =  (N-l)(N-2)/2  >  0  as  long  as  N  >  2,  so  it  would  be  applicable  for  smaller  N 
provided  the  assumption  (2.24)  can  be  made  and  predictions  beyond  j  =  N  are  not  needed. 


Continuing  with  the  development  of  a  prediction  interval,  (2.33)  leads  to 


p|exp(pjjtj^2^"-(N-l)))  <  <  exp(pj.tj  ^2(''-(N-1)))  }  =  1-a,  (2.34) 

which  in  turn  implies  that 

p{c*jexp(pjj5  lo^2(n-(N-l)))<  Cjj  <  C*.exp(pj.5  }  =  l-o.  (2.35) 

Relation  (2.35)  says  that  with  probability  1-a,  the  future  cost  C. .  will  lie  in  the  interval 

^«J 

/c*jexp(pjja  1))) .  C*jexp(pjja  (2.36a) 


In  (2.36a),  it  is  tacitly  assumed  that  0  <  a  <  1.  Notice  that  the  left  endpoint  is  the  point 

A  ♦ 

prediction  C-  j  multiplied  by  a  factor  less  than  1  (i.e.  reduced,  since  tim)  <  0  for  u  <  1/2)  and 

*>J  D 

A  ♦ 

the  right  endpoint  is  the  point  prediction  C. .  multiplied  by  a  factor  greater  than  1  (i.e.  increased, 
since  t.^(m)  >  0  for  u  >  1/2). 

By  taking  C-  ,  s  1  for  all  i,  the  analysis  that  led  to  expression  (2.36a)  also  provides  (1-a) 
prediction  intervals  for  cumulative  LDF  via 


(“P&k=2^ik  +  Pij®  •  “p(4=2^ik  Pij*l-a«("-(N-'») 

Similarly,  prediction  limits  for  individual  LDFs  are  easily  developed.  Letting 


(2.36b) 


s?  =  Wj  +  Wj/(N-j+l) 


(2.36c) 
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it  follows  that  a  1-a  prediction  interval  for  the  LDF  exp(Yjj)  is  (assuming  j  >  N-i+1) 

|exp(Y*.  +  s.a  .  exp(Yjj  +  s-o  t, 

As  before,  the  prediction  formula  can  be  expressed  as 


(2.36d) 


A  *  A  *  A  * 


When  N-i+1  <  j  <  N,  (2.25)  and  (2.27a)  imply  that 


(2.37) 


A 

C 


1,J  i.J-l\llr=l  rj; 


(2.38) 


where 


LDFjj  =  exp(Yj.j)  (2.39) 

is  the  observed  loss  development  factor  from  year  j-1  to  j  for  cohort  r.  Thus,  this  method  of 
prediction  is  nearly  the  same  as  the  basic  actuarial  method,  except  that  simple  averaging  of  LDFs 
has  been  replaced  by  geometric  averaging  of  the  LDFs,  as  seen  in  (2.38). 
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3.  Application  to  U.S.  Department  of  the  Navy  Workers’  Compensation  Claims 

3.1  Application  to  Complete  1990  Through  1993  Cohort  Data 

The  data  in  Table  1.1  of  the  introduction  represent  actual  complete  cohort  data  for  1990 
through  1993  cohorts  having  workers’  compensation  claims  against  the  U.S.  Department  of  the 
Navy.  This  data  set  is  small,  and  the  LDF  trend,  if  any,  across  cohorts  is  weak,  as  seen  in  Table 
1.2.  The  full  parametric  model  of  sections  2.2  through  2.3  is  not  applicable  (too  many  unknown 
parameters  in  that  model),  but  the  model  of  section  2.4  may  be  applied  to  compute  predictions 
and  prediction  intervals.  Using  N  =  4,  w^  =  ...  =  w^  =  1,  formulae  (2.28)  and  (2.31)  yield  a  = 

0.0299184;  P24  =  2,  P33  =  3/2,  P34  =  7/2,  =  ^^3,  P43  =  17/6,  and  P44  =  29/6.  Here, 

(N-l)(N-2)/2  =  3.  Taking  a  =  .20  for  an  80%  prediction  interval,  t  9q(3)  =  1.6377  =  -t  jq(3),  and 

(2.36)  yields  the  prediction  intervals  in  Table  3.1. 


Table  3.1.  80%  Prediction  Intervals  for  1991  Through  1993  Cohorts  Based  on  Data  from 


Table  1.1  (In  Thousands  of  Dollars) 

Develooment  Year 

Cohort  Year 

2 

3 

4 

Lower 

Pred. 

UPBgr 

Lower  Pred.  Unoer 

Lower 

Pred.  Uoner 

1991 

71,343 

76,462  81,948 

1992 

59,375  63,047  66,946 

73,140 

80,161  87,857 

1993 

37,263 

39,432  41,727 

54,646  59,344  64,446 

67,749 

75,454  84,036 

3 

3.2  Application  to  Actuarial  Data  From  Miccolis 

The  somewhat  wide  prediction  intervals  in  Table  3.1  reflect  the  fact  that  they  are  based 
on  a  very  small  sample.  When  more  data  arc  available,  the  regression  model  of  sections  2.1 
through  2.3  can  lead  to  more  precise  predictions.  To  illustrate  their  use,  these  models  were 

-j 

applied  to  data  found  in  exhibit  8,  sheets  3a  and  3b  of  the  report  by  Miccolis.  These  data 
represent  (mostly  imputed)  total  cumulative  costs  (indemnity  plus  medical)  for  cohorts  from 
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1961  through  1993,  and  development  years  1  to  32.  Because  complete  data  were  unavailable  for 

3 

cohort  years  prior  to  1990,  Miccolis  employed  a  variety  of  actuarial  techniques  to  estimate 
claim  counts  and  average  costs  per  claim  for  the  cohorts  prior  to  1990  using  all  available  data. 
Thus,  the  reader  shou''^  keep  in  mind  that  data  employed  in  this  example  contain  a  great  deal  of 
imputed  values  and  are  smoother  (i.e.  exhibit  less  fluctuation)  than  a  completely  real  data  set 
Consequently,  the  variance  estimate  (2.10)  will  be  smaller  than  expected,  leading  to  prediction 
intervals  that  are  narrower  than  could  be  expected  from  comparable  real  data. 

For  this  application,  there  are  N  =  33  cohorts,  and  the  cumulative  cost  data  are 
summarized  in  Table  3.2.  The  data  employed  by  the  models  of  sections  2. 1  through  2.3  are  the 
log-LDFs,  namely  Y--  =  log(C-  -/C-  •  ,),  1  <  i  <  34-j,  2  <  j  <  32,  C-  •  being  the  cumulative  cost 

of  cohort  i  through  development  year  j.  Cohort  years  1961  through  1993  correspond  to  i  =  1 
through  i  =  33,  respectively. 
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Table  3.2.  Cumulative  Cost  Data  From  Miccolls 
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Table  3.2,  continued.  Cumulative  Cost  Data  From  Miccolis 


Using  stepwise  regression,  the  model  (2.6)  was  fit  without  weighting  (i.e.  using  Wj^=l  for 
all  k).  Only  seven  terms  in  (2.5)  were  needed,  and  the  best  model  therefore  had  p  =  7  and 

Yij  =  v(ij)‘p  +  e.j  (3.1) 

with 

v(i,j)‘  =  (l  I/j^  i/j^  1/)“  i/j'*  iV  (3.2) 

and 

P^  =  (Aoo  Aq2  Aq3  Aq^  a A22  A23).  (3.3) 

2 

The  estimates  of  P  and  a  were 

^^  =  (.0214289  8.6704  -22.5927  31.5725  -.016130  -00123026  .0005548) 

(3.4) 

=  (.007093)^ 

and  the  fit  yielded  an  value  of  99.9%. 

Miccolis  provided  predictions  of  the  costs  C-  •  for  35-i  <  j  <  32,  3  <  i  <  33.  Converting 
those  and  the  data  C. .,  1  <  i  <  34-j,  1  <  j  <  32,  by  the  formula  LDF..  =  C-  /C-  •  , ,  and  plotting 
LDFy  versus  (i,  j)  yields  Figure  3.1.  Using  the  fitted  model  (3.1)  through  (3.4),  and  using  the 

A 

methods  in  sections  2.2  to  compute  predictions  C-  •,  the  surface  in  Figure  3.2  was  produced.  It  is 

readily  seen  that  not  only  does  the  model  (3.1)-(3.4)  provide  an  excellent  fit  to  the  data,  but  it 
also  produces  predictions  of  the  LDFs  (outside  of  the  original  data)  that  are  close  to  the 
predictions  produced  by  Miccolis  using  different  actuarial  methods.  It  is  emphasized  here  that 
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3 

the  predicted  LDFs  from  Miccolis  (lower  right  part  of  the  data  matrix  of  Table  3.2)  were  NOT 
part  of  the  data  used  to  fit  the  regression  model  (3.1)-(3.4). 
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Figure  3.1.  Loss  Development  Factors 
(Observed  and  Predicted  From  Miccolis3) 


Using  (2.23b)  and  (2.36b),  predictions  and  prediction  interval  endpoints  were  computed 
for  cumulative  LDF  (CLDF)  and  LDF.  (Recall  that  an  LDF  is  always  for  consecutive  years  in 
this  report,  while  cumulative  LDF  describes  the  loss  development  from  the  first  year  to  the  year 
specified.)  These  are  shown  in  Tables  3.3  and  3.4.  Table  3.5  shows  all  the  predicted  costs  up  to 
year  32  for  the  1989  through  1993  cohorts,  along  with  lower  and  upper  endpoints  for  prediction 
intervals.  Figure  3.3  shows  the  cumulative  cost  predictions  for  the  1990  cohort,  along  with  upper 
and  lower  endpoints  of  a  95%  prediction  interval.  The  costs  prior  to  year  5  are  known  for  the 
1990  cohort,  so  there  are  no  prediction  intervals  prior  to  year  5. 

Recall  that  the  LDF  and  CLDF  are  known  exactly  up  to  and  including  years  5,  4,  3,  and  2 
for  the  cohorts  of  years  1989,  1990,  1991,  and  1992,  respectively.  Thus,  these  are  not  tabled. 
Also,  the  cumulative  LDF  at  year  5,  4,  3,  and  2  has  been  divided  out  of  the  formulae  (2.23b)  and 
(2.36b)  for  cohorts  1989,  1990,  1991,  and  1992,  respectively.  That  is,  the  CLDFs  in  Table  3.3 
start  with  the  initial  year  shown  as  the  first  year  in  the  cumulative  calculation.  Thus,  to  use  the 
tables  to  compute,  for  example,  the  predicted  cost  at  development  year  32  for  the  1989  cohort, 
look  up  the  CLDF  for  year  32  in  Table  3.3  and  multiply  it  by  the  last  known  cumulative  cost, 
namely,  that  at  year  5.  This  value  is  (in  thousands)  $86,373  (from  Table  3.2).  Thus,  the 
prediction  is  (4.2077)  ($86,373)  =  $363,432  which  is  in  agreement  (within  round-off  error)  with 
Table  3.5.  This  usage  of  the  term  cumulative  LDF  is  slightly  different  from  common  usage  in 
actuarial  literature.  See  section  2.2  for  a  discussion  of  these  differences. 

As  another  example,  to  compute  a  prediction  and  prediction  interval  for  the  cumulative 
cost  for  the  1993  cohort  at  year  32,  look  up  the  predicted  CLDF  from  Table  3.3,  Cohort  Year  = 
1993,  and  Year  =  32.  The  predicted  value  is  26.8176.  For  a  95%  prediction  interval  (with  equal 
tail  probabilities),  find  the  .025  and  .975  endpoints  for  year  32.  These  are  24.6378  and  29.1903, 
respectively.  Finally,  multiply  these  (the  endpoints  and  the  prediction)  by  the  initial  (year  1) 
cost  for  the  1993  cohort  (from  Table  3.2,  $13,702  in  thousands).  For  this  example,  the  predicted 
cost  year  32  is  then  $367,455,  and  the  95%  prediction  interval  is  ($337,587,  $399,965),  all  in 
thousands. 
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Figure  3.3.  95%  Prediction  Intervals 
1990  Cohort  (Costs  in  $1 K) 


Table  3J.  Cumidative  Loss  DevchqHncnt  Factors:  Predicted  and  Prediction  Interval  Endpoints 


Table  Cumulative  Lon  Dcvdopment  Factors:  Predicted  and  Prediction  Interval  Endpoints 


Tabic  3w3.  Cumulative  Loss  Development  Factors:  Predicted  and  Prediction  Interval  Endpoints 
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*  These  apply  to  the  1993  cohort  only.  While  numbers  for  cohort  years  in  the  near  future 
will  be  similar,  caution  should  be  used  in  basing  predictions  for  cohorts  in  the  distant  future 
on  these  for  the  1993  cohort. 

**  This  can  be  interpreted  as  a  $1 .00  cost  in  year  1 .  For  example,  an  initial  cost  of  $1 .00 
in  1993  will  grow  to  a  total  accumulated  cost  of  $26.82  by  year  32. 


Table  3.4.  Loss  Development  Factor:  Predicted  and  Prediction  Interval  Endpoints 


Table  3.4.  Loss  Developmeot  Factor:  Predicted  and  Prediction  Interval  Endpoints 
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TabI  j  3^.  Loss  Development  Factor:  Predicted  and  Prediction  Interval  Endpoints 


1  Coh«i¥i^t9»  1 

0.0100 

O.<B50 

0.0500 

0.9000 

Q.9900 

0.9750 

2.7918 

2.8001 

2.8073 

2.8155 

2.8450 

2.8747 

2.8831 

2.8905 

2.8991 

1.4728 

1.4769 

1.4804 

1.4845 

1.4990 

1.5135 

1  5177 

1.5213 

1.5255 

1.2714 

1.2748 

1.2778 

1.2812 

1.2934 

1.3057 

1.3092 

1.3122 

1.3158 

r  * 

1.1871 

1.1902 

1.1930 

1.1961 

1.2073 

1.2187 

1.2219 

1.2247 

1.2280 

i  6 

1.1393 

1.1423 

1.1449 

1.1479 

1.1586 

1.1694 

1.1724 

1.1751 

1.1782 

. t 

1.1087 

1.11)6 

1.1141 

1.1170 

1.1273 

1.1378 

1.1407 

1.1433 

1.1463 

_ 

1.0876 

1.0904 

1.0929 

1.0957 

1.1058 

1.1160 

1.1189 

1.1214 

1.1244 

1.0724 

1.0752 

1.0776 

1.0804 

1.0904 

1.1004 

1.1032 

1.1057 

1.1086 

— ifn 

1.0611 

1.0638 

1.0662 

1.0690 

1.0788 

1.0887 

1.0915 

1.0940 

1.0968 

1.0524 

1.0551 

1.0575 

1.0602 

1.0700 

1.0798 

1.0826 

1.0850 

1.0878 

M^KI 

1.0483 

1.0506 

1.0534 

1.0630 

1.0728 

1.0755 

1.0779 

1.0808 

1.0401 

1.0428 

1.0452 

1.0479 

1.0575 

1.0672 

1.0699 

1.0723 

1.0751 

m 

1.0357 

1.0384 

1.0407 

1.0434 

1.0530 

1.0626 

1.0654 

1.0678 

— iq 

1.0321 

1.0348 

1.0371 

1.0398 

1.0493 

1.0589 

1.0616 

1.0640 

mmti 

1.0290 

1.0317 

1.0340 

1.0367 

1.0462 

1.0558 

1.0585 

1.0609 

1.0637 

[ _ II 

1.0265 

1.0292 

1.0315 

1.0341 

1.0436 

1.0532 

1.0559 

1.0583 

1.0610 

m 

liWiri 

1.0243 

1.0270 

1.0293 

1.0320 

1.0414 

1.0509 

1.0537 

1.0560 

1.0588 

1.0225 

1.0251 

1.0274 

1.0301 

1.0395 

1.0490 

1.0518 

1.0541 

1.0569 

1 _ w 

1.0209 

1.0235 

1.0258 

1.0285 

1.0379 

1.0474 

1.0501 

1.0525 

1.0552 

1.0195 

1.0221 

1.0244 

1.0271 

1.0365 

1.0460 

1.0487 

1.0510 

1.0538 

1.0183 

1.0209 

1.0232 

1.0259 

1.0352 

1.0447 

1.0474 

1.0498 

1.0525 

1.0172 

1.0198 

1.0221 

1.0248 

1.0342 

1.0436 

1.0463 

1.0487 

1.0514 

[  . ^ 

1.0162 

1.0189 

1.0212 

1.0238 

1.0332 

1.0427 

1.0454 

1.0477 

1.0504 

1.0154 

1.0181 

1.0203 

1.0230 

1.0323 

1.0418 

1.0445 

1.0468 

1.0496 

^^173 

1.0147 

1.0173 

1.0196 

!.0222 

1.0316 

1.0410 

1.0437 

1.0461 

1.0488 

iHm 

1.0140 

1.0166 

1.0189 

1.0215 

1.0309 

1.0403 

1.0430 

1.0454 

1.0481 

1.0134 

1.0160 

1.U183 

1.0209 

1.0303 

1.0397 

1.0424 

1.0447 

1.0475 

1 _ » 

1.0128 

1.0155 

1.0177 

1.0204 

1.0297 

1.0392 

1.0418 

1.0442 

1.0469 

MiKIO 

1.0123 

1.0150 

1.0172 

1.0199 

1.0292 

1.0386 

1.0413 

1.0437 

1.0464 

1.0119 

1.0145 

1.0168 

1.0194 

1.0288 

1.0382 

1.0409 

1.0432 

1.0459 

rnmmm^ 

1.0115 

1.0141 

1.0164 

1.0190 

1.0283 

1.0378 

1.0404 

1.0428 

1.0455 
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Table  3.5.  Cumulative  Costs:  Predicted  and  Prediction  Intervai  Endpoints  (In  Thousands  ot 


0.0100 


99.196 


111.691 


123.394 


134.404 


144.829 


154.767 


164.300 


173.502 


182.432 


191.141 


199.673 


208.065 


216.349 


224.554 


232.703 


240.818 


248.918 


257.018 


265.135 


273.282 


281.471 


289.713 


298.018 


306.397 


314.857 


323.408 


332.057 


99.456 


112.109 


123.963 


135.122 


145.696 


155.783 


165.467 


174.821 


183.903 


192.767 


201.456 


210.008 


218.454 


226.823 


235.140 


243.425 


251.699 


259.977 


268.276 


276.608 


284,987 


293.423 


301.927 


310.510 


319.179 


327,943 


336.812 


0<3500 


99.681 


112,470 


124.454 


135,742 


146.446 


156.663 


166.478 


175,963 


185,179 


194.178 


203.003 


211.693 


220.281 


228.793 


237.256 


245.691 


254,116 


262.550 


271.007 


279.502 


288.046 


296,652 


305.331 


314.092 


322.944 


331.896 


340.956 


125.023 


136.462 


147.315 


195.817 


204.802 


231.086 


mmsmi 


248.329 


265.548 


282.875 


291.615 


300.420 


309.303 


318.273 


■ESESII 


100.864 


114.374 


127,052 


139.029 


150.423 


161.336 


171.853 


182.048 


191.982 


201.709 


211,274 


220,716 


230.068 


239.360 


248.617 


257.863 


267.118 


276.399 


285.723 


295.106 


304,560 


314,098 


323.732 


333.473 


343,331 


353.315 


363.435 


0.9000 


101.794 


115.880 


129.113 


141.645 


153.597 


165.074 


176,161 


186.934 


197,456 


207.779 


217.951 


228.011 


237.994 


247.930 


257.845 


267,764 


277.707 


287.694 


297.741 


307.865 


318.080 


328.399 


338.834 


349,399 


360.103 


370,957 


381.971 


0  9500 


102,060 


116.311 


129,704 


142.395 


154.509 


166.149 


177.402 


188.343 


199.036 


209.533 


219.882 


230.123 


240.290 


250.414 


260.523 


270.639 


280.784 


290.979 


301.239 


311.581 


322.020 


332.570 


343.242 


354,050 


365.004 


376.116 


387.395 


0.9730 


102.291 


116.686 


130.218 


143.049 


155.304 


167.087 


178.486 


189,574 


200.416 


211.066 


221.570 


231.970 


242.299 


252.589 


262.867 


273,157 


283.481 


304.306 


314.841 


325.477 


336.230 


347.112 


358.135 


369.310 


380.649 


392.162 


0:9900 


102.560 


117.123 


130.818 


143.813 


156.234 


168.184 


191.015 


202.033 


212.862 


223.550 


234.136 


244.656 


255.142 


265.620 


276.115 


286.649 


307,910 


318.673 


329,543 


340,536 


351.665 


mUMiMX 


0.0t0l} 


95,758 


110.811 


124.778 


137.807 


150.045 


161.626 


172,660 


183.246 


193.461 


203.376 


213.046 


222.520 


231.840 


241.041 


!■E3a]^l 


259.209 


268.225 


277.225 


286.228 


295.249 


304.304 


313.407 


322.570 


331.805 


341.121 


350.529 


360.038 


369.656 


0.0250 


96.011 


111.230 


125.361 


138.554 


150,957 


162,704 


173.906 


184,659 


195.045 


205.130 


214.974 


224.624 


234,122 


243.504 


252.802 


262.043 


271.250 


280.445 


289.646 


298.871 


308.133 


317.448 


326.828 


336,285 


345.829 


355.470 


365.218 


375.081 


96,230 


111.593 


125.865 


139,200 


151,746 


163.637 


174,984 


185.884 


196.417 


206.652 


216.646 


226.449 


236.103 


245,643 


255.102 


264.506 


273.880 


283,245 


292.619 


302.021 


311.466 


320,966 


330.536 


340.187 


349.931 


359.776 


369.733 


379.812 


96,482 


112.012 


126.449 


139.948 


152.661 


164.719 


176.235 


187.306 


198.011 


208.420 


218.591 


228.573 


238.408 


248.133 


257.779 


267.374 


276.943 


286,507 


296.086 


305.695 


315.352 


325.071 


334.864 


344,743 


354.720 


364.806 


375.009 


385.340 


y6*ft990 


Actual 


14.955 


41.424 


62.897 


79.971 


Pliiilicied 


97.379 


113.503 


128.528 


142.621 


155.933 


168.595 


180.722 


192.410 


203.739 


214.780 


225,592 


236.224 


246.720 


257,118 


267.450 


277.744 


288.027 


298,320 


308.644 


319.017 


329.455 


339,973 


350.586 


361.307 


372.146 


383.117 


394.229 


405.492 


09000 


98.284 


115.015 


130.642 


145.345 


159,275 


172.563 


185.323 


197,653 


209.633 


221.334 


232.816 


244.131 


255,322 


266.428 


277.483 


288.516 


299.554 


310.620 


321.735 


332.919 


344.188 


355.559 


367.047 


378.666 


390,429 


402.347 


414,433 


426.698 


98,542 


115.447 


131.248 


146.126 


160.235 


173.704 


186.649 


199,165 


211,334 


223.228 


234.906 


246.420 


257,814 


269.128 


280.395 


291.645 


302,905 


314.198 


325.546 


336.968 


348.483 


360.106 


371.853 


383,737 


395.773 


407,972 


420,347 


432.909 


O.97S0 


98,766 


115,823 


131.775 


146.808 


161,072 


174,701 


187.806 


200.485 


212.821 


224.884 


236.734 


248.423 


271,492 


282.945 


294,386 


305.841 


317,335 


328.888 


340,521 


352.252 


364,097 


376.071 


388.190 


400,467 


412.914 


425.544 


438.368 


0.9900 


99.028 


116.261 


132.391 


147.604 


162,051 


175.866 


202.032 


214.563 


226.824 


238.876 


250.772 


262.555 


274,266 


285.939 


297.605 


309.291 


321.020 


332.816 


344.698 


356.684 


393.432 


405.994 


418.734 


431.667 


444.802 
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Table  3.5.  Cumulative  Costs: 


Predicted  and  Prediction  Interval  Endpoints  (In  Thousands  of  $) 


CohonY«tf)9»l 


(LOtOO 


77,151 


92,990 


107.551 


121,018 


133,566 


145.345 


156,489 


167,107 


177,293 


187.124 


196,665 


205.973 


215.094 


224.068 


232.930 


241,708 


250.429 


259,116 


267.788 


276.463 


77,358 


93,349 


108.065 


121.690 


134.398 


146.339 


168,427 


178.778 


188,775 


198.484 


MKMSMI 


293.887 


302,662 


311,496 


IKSIEEII 


329,383 


338.455 


347.625 


356.902 


217.254 


226.402 


235,439 


244.397 


253.300 


262.173 


271.034 


279.903 


288,795 


297,725 


306.707 


77.536 


93 


108,509 


122.271 


135.118 


147,200 


158,647 


169.571 


180.065 


190.207 


200.062 


209.687 


219,129 


228.428 


237.620 


246.734 


255.796 


264.831 


273,858 


282.896 


170.900 


\wmnmMi 


201.897 


211.694 


286.387 


II 


324.871 


334.075 


343.374 


352.777 


362.292 


II 


301.068 


310,230 


319.459 


328.768 


338.166 


347.663 


357.270 


366.993 


314.341 


■iSSEOII 


342.944 


352.675 


362.521 


372.489 


95.292 


110.856 


125.350 


138.940 


151.774 


163.981 


175.671 


186.938 


197.860 


208.505 


218.930 


229.183 


239.306 


249.334 


259.301 


269.233 


279.154 


289.086 


299.048 


309,059 


319.134 


329.287 


339.532 


349.881 


360.347 


370,939 


381.668 


392.544 


q;9Qoo 


79.211 


96.585 


112.720 


127.802 


141.992 


155.436 


168.260 


180.575 


192.474 


204.036 


215.330 


226.413 


237.335 


248.138 


258.861 


269.535 


280.188 


290.847 


301.535 


312.270 


323,072 


333.958 


344.943 


356.042 


367.267 


378,632 


390.148 


401.827 


413.678 


a9500 


79.421 


96.954 


113.254 


128.506 


142.870 


156,490 


169.494 


181.990 


194.073 


205.822 


217.305 


228.580 


239.698 


i^^EQll 


261.627 


272.508 


283.374 


294.251 


305.160 


316.123 


327.159 


338.284 


349.515 


360.866 


372,351 


383.982 


419.874 


a9750 


79.604 


97.276 


113,719 


129.119 


143.635 


157,410 


170,571 


183.226 


195.470 


207.383 


219.032 


230.477 


241.766 


252.945 


264.050 


275.114 


286.167 


297,234 


308.340 


319.504 


342.082 


353,529 


365.103 


376.817 


388.683 


400.716 


412.925 


425.323 


62.189 


80.168 


96.485 


111.448 


125,275 


138.151 


161.669 


172.562 


183.013 


193.100 


202.892 


212.446 


221.809 


240.121 


249.136 


258.094 


267.018 


275.927 


amso 


62,360 


80,487 


96.962 


1 12.090 


126.084 


139.130 


151,387 


162.990 


174.056 


184.679 


194.941 


204.909 


214.640 


233.578 


242.861 


252,064 


261.212 


279,437 


288.553 


320.848 


330.001 


339.237 


348.565 


357.995 


367.535 


316.110 


325.410 


334.788 


344,255 


353.819 


363.490 


373.278 


62,507 


80.763 


97,375 


112.645 


126.784 


139,977 


152.383 


164.136 


175.351 


186.125 


196.539 


206.660 


216.546 


226.245 


235,799 


245,243 


254.609 


263.924 


273,211 


282.492 


291,784 


301.105 


310.470 


319.893 


329.387 


338.962 


348.630 


358.401 


368.285 


378.290 


97.853 


113.288 


127.597 


140,961 


153,539 


176,856 


187.807 


198.398 


208.698 


218.764 


228.646 


238.386 


248.019 


267.086 


276,572 


286.055 


295.554 


305,087 


314.668 


334.031 


363,758 


■KSHnil 


y«arl992 


AfiiuU 


14.468 


41.982 


63.282 


82.219 


99.557 


115,585 


130.503 


144.486 


157.689 


170.246 


182.270 


193.859 


205.095 


216.047 


226.773 


237.323 


247,741 


258,063 


268.322 


278.546 


288.761 


298.988 


309.247 


319,557 


329,934 


340,392 


350.946 


361.608 


372.390 


383.304 


394.360 


405,567 


&9O0O 


63.893 


83,371 


101.291 


117.930 


133.476 


148.100 


161.951 


175,163 


187.850 


200.108 


212.019 


223,655 


235.074 


246.329 


0.9500 


64,067 


118.603 


134.331 


149.141 


268.514 


279.516 


290.498 


301.487 


312.506 


323.575 


334.714 


345.941 


357.270 


368.717 


380.296 


392.020 


403,900 


415.948 


428,176 


225.860 


237,483 


248,944 


260.287 


271.552 


293.978 


305.195 


316.447 


350.618 


362.205 


385,767 


397,770 


422.280 


434.811 


64.219 


83,987 


102.221 


119.190 


135.077 


150.049 


164.254 


177.824 


190.872 


203.496 


215.778 


227,790 


239.591 


251.234 


262,762 


274.216 


285.629 


297.031 


308.449 


319.907 


331.426 


343.026 


354.726 


366.540 


378.486 


390,576 


402.826 


415.246 


427.850 


440.650 


79.817 


97.651 


114.263 


129.836 


144.530 


158.487 


171.831 


184.674 


197,108 


209.213 


221.058 


232 


244.194 


255.579 


266.895 


278.175 


289.448 


300.741 


312.078 


323.479 


EEnl 


346.550 


358.254 


370.091 


382.076 


394.221 


406.540 


419.044 


431.745 


64.395 


84.322 


102.727 


119.876 


135.950 


151,112 


165.511 


179.277 


192.524 


205.349 


217.836 


230.055 


242.066 


253.922 


265.669 


277,345 


288.985 


300.619 


312.274 


323,975 


335.744 


347.600 


359.561 


371.645 


383.867 


396,242 


408.784 


421.505 


434.418 


447.535 


Table  3.5.  Cumulative  Costs:  Predicted  and  Prediction  Interval  Endpoints  (In  Thousands  ofS) 


38.253 


56.912 


73.178 


87.907 


101.400 


113.860 


125.460 


136.345 


146.641 


156.452 


165.865 


174,952 


183.773 


192.381 


200.819 


209.123 


217,325 


225.453 


233.530 


241,577 


249.613 


257,654 


265.714 


273.808 


281.947 


290.141 


298.401 


306.737 


315,157 


323,669 


332,280 


00230 


38.366 


57.148 


73,550 


88.425 


102.068 


114.681 


126.436 


00300 


38.465 


57.353 


73.872 


88.872 


WBSSSMl 


157,898 


167.469 


176.714 


185.696 


194.465 


203.066 


211.535 


219.905 


228.203 


236,453 


244.676 


252.892 


261.116 


269.363 


277.648 


285.982 


294.376 


302.840 


3U.385 


320.019 


328.750 


337.587 


149.048 


196.276 


205.019 


213.633 


230,595 


238,997 


247.374 


255.747 


264.131 


272.542 


280.994 


289,499 


298.068 


306,711 


315.440 


324.262 


333.185 


342,219 


38.578 


57.589 


74.245 


89.391 


103.316 


116.218 


128.264 


139.599 


150,347 


160.612 


170.480 


180.025 


189.310 


198.385 


207.294 


216,077 


233.385 


241.964 


250.522 


259.079 


284.902 


293.608 


302.382 


311.236 


320,180 


329.223 


338.373 


347.640 


38.981 


58.431 


75.575 


91.244 


105.716 


119.177 


131.790 


143.698 


155.021 


165.865 


176.317 


186.451 


196.329 


206,005 


215.524 


224.924 


252,727 


261.949 


271.182 


280.446 


289.756 


299.128 


308.574 


327.738 


337,479 


347.339 


357.328 


367.455 


39.388 


59.286 


76,929 


93.137 


108.171 


122.211 


135.413 


147.916 


159.841 


171.291 


182.354 


193.106 


203.610 


213.919 


224.080 


234.133 


263.970 


273,897 


283.851 


293,853 


0.9500 


39,504 


59,530 


77,318 


93.680 


108.877 


123.085 


136,459 


149,135 


345.115 


355.712 


366.452 


388.399 


172.861 


184.103 


195,036 


205.722 


216,217 


226.567 


236.812 


246.987 


257.123 


267.247 


277.382 


287.549 


297.769 


308.058 


318.432 


328.905 


339.493 


350.206 


361.058 


372.058 


383.219 


394.550 


0.9750 


39.605 


59.743 


77.656 


94.154 


109.494 


123.849 


137.371 


150.199 


162,452 


174.234 


185.633 


196.725 


207.572 


218.230 


228,746 


239.160 


249,507 


270.122 


280.440 


290.796 


301.207 


311.693 


322.269 


332.950 


343.751 


354,683 


365.759 


376.991 


388.389 


09900 


39.723 


59.991 


78,051 


94.708 


110.215 


124.742 


138.440 


151.447 


163.880 


175.845 


187.428 


198.706 


209,743 


220.594 


231.305 


241.919 


252.469 


262.988 


273.503 


284.038 


294.615 


305.254 


315.973 


326,789 


337,715 


348.768 


359.959 


37I.30I 


382.807 


394.487 


406.352 
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4.  Conclusions  and  Recommendations 


This  investigation  shows  that  the  Loss  Development  Method  commonly  used  by  actuaries 
to  predict  workers’  compensation  costs  can  be  cast  in  the  context  of  intrinsically  linear  models 
and  thereby  made  amenable  to  the  theory  of  linear  statistical  models  for  the  computation  of  point 
and  interval  predictions  of  future  costs.  These  computations  have  been  illustrated  using  actual 
and  imputed  U.S.  Department  of  the  Navy  workers’  compensation  claims.  In  addition,  the 
log-loss  development  regression  model  developed  herein  has  been  shown  to  produce  point 
predictions  that  are  nearly  the  same  as  those  produced  (with  apparently  much  more  effort)  using 
traditional  actuarial  methods  (Figures  3.1  and  3.2).  In  contrast  to  the  actuarial  approach,  the 
regression  approach  is  relatively  easy  to  compute,  with  the  computations  involved  being 
standard  in  many  statistical  computer  packages,  and  it  provides  a  means  of  assessing  the 
accuracy  of  the  resultant  predictions.  Moreover,  beyond  knowledge  of  basic  linear  statistical 
models  and  statistical  analysis,  specialized  knowledge  is  not  needed  to  apply  them. 

I  emphasize  here  that  the  reader  should  not  conclude  that  an  inexperienced  analyst  armed 
with  the  linear  model  methods  employed  herein  can  somehow  replace  the  traditional  analyses 
and/or  services  of  a  professional  actuary.  Rather,  the  methods  of  this  investigation  should  be 
evaluated  further  by  actuarial  scientists  and  practitioners  and  perhaps  be  adopted  (with  any 
necessary  and  suitable  modifications)  as  another  set  of  tools  in  the  collection  of  numerical 
methods  that  have  come  into  actuarial  practice. 

Mathematically  and  statistically  speaking,  a  few  paths  for  future  research  should  be 
pursued.  First,  recall  that  the  linear  model  formulated  in  section  2  assumes  that  the  error 
variances  are  not  all  equal.  In  fact,  they  are  assumed  to  approach  zero  as  development  year 
increases.  This  is  a  reasonable  assumption  and  one  that  should  be  retained.  However,  the 
examples  used  to  illustrate  the  methods  all  assumed  homogeneity  of  variances.  This  approach 
was  taken  in  section  3.1  because  the  data  set  was  so  small  and  in  section  3.2  because  the  data 
(which  contained  a  great  deal  of  imputed  data)  was  already  very  smooth  and  the  errors  very 
small.  In  practice,  an  effort  should  be  made  to  assess  the  weights  in  (2.2)  perhaps  considering 
expansions  of  Wj  that  would  lead  to  Wj  =  0(l/j'^)  as  for  some  v  >  1. 
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A  crucial  assumption  that  allows  the  development  of  closed-form  expressions  for  the 
prediction  interval  endpoints  is  that  of  a  Gaussian  error  distribution.  Some  investigation  should 
be  undertaken  using  more  extensive  data  to  test  the  Gaussian  assumption.  This  was  not  pursued 
in  this  investigation  because  the  example  of  section  3.1  had  too  few  data  points  and  the  data 
from  section  3.2  had  too  much  imputed  data. 

In  parallel  with  testing  the  Gaussian  assumption,  alternative  methods  (not  based  on  the 
Gaussian  assumption)  of  computing  prediction  intervals  should  be  investigated.  In  particular, 
the  models  of  section  2  are  all  amenable  to  bootstrapping,  a  computational  method  that  allows 
the  estimation  of  the  relevant  statistical  quantities  with  only  minimal  assumptions.  I  am 
currently  pursuing  this  in  a  separate  investigation  that  will  make  use  only  of  the  actual  U.S. 

3 

Department  of  the  Navy  data  (as  opposed  to  the  imputed  data  in  the  study  by  Miccolis  ). 

Finally,  I  recommend  that  models  having  different  regression  structures  (perhaps 
nonlinear,  and  not  intrinsically  linear),  as  well  as  some  classes  of  dynamic  statistical  models  (for 
example  e^),  t  =  1,  2,  ...  )  be  investigated.  The  actuarial  literature  is  currently 

focused  on  ad  hoc  (although  effective  and  useful)  curve  graduation  and  point  projection  methods 
and  is  relatively  short  of  models  with  sufficient  structure  to  allow  probabilistic  assessment  of 
prediction  accuracy.  This  report  provides  a  stimulus  in  these  directions. 
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